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Base on the concept of complementary media, a scheme to achieve cloak an object outside the 
shell is proposed by Yun lai et al [PRL 102, 093901] .We generalize the concept in acoustic waves 
and demonstrate numerically how to achieve the acoustic cloaking and imaging with complementary 
media. 

PACS numbers: 43.20.Fn 

Making use of form-invariant coordinate transformation of Maxwell's equation, J. B. Pendry and Leonhardt inde- 
pendently proposed the invisibility cloak and the concept of transformation media for electromagnetic (EM) waves 
[U [2]. Transformation media transport electromagnetic (EM) waves without reflection, giving rise to exotic EM 
devices such as wave concentrator's , cylindrical superlens|4|, etc., besides invisibility cloak [5l El III [HI 1^. From the 
view of transformation optics, left-handed materials are precisely the complementary media of conventional dielectric 
materials by mirror transformation pUl llll [T^ I13j . Under appropriate design, a complementary object destructively 
interferes with an outside object, and perfectly cancels the scattered waves from it, giving rise to perfect cloak to the 
whole system[T4]. The idea is also successfully applied to design perfect lens [15], the "superscatterer" fT6| and the 
"anticloak" [T7J (TH] as well. The concept of transformation media can be extended to the scope of acoustic waves. 
A scheme to achieve acoustic cloaking in two dimensions was proposed by Cummer and SchurigfTO]. Subsequently, 
Huanyang Chen and C. T. Chan generalized the idea to three dimensions by mapping the acoustic equation directly 
to the conductivity equation [20J. Recently, a technique to achieve broadband acoustic cloak was proposed by J. B. 
Pendry and Jensen Li|21j. 

In this paper, we show that the concept of complementary media can be generalized in acoustic waves. The 
acoustic complementary media possess constitute parameters in opposite sigh of each other. Perfect acoustic cloaking 
and perfect imaging with complementary media are demonstrated numerically in two dimensions. The recovery of 
the acoustic wave field by complementary media is interpreted in one dimensional multilayer system. 

The acoustic wave equation takes a form as: 
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where mass density p(a;)and bulk modulus k{x) are two constitutive parameters of the acoustic media at x point. 
An alternative parameter equivalent to k {x) is sound velocity c (x) with {x) = k (x) / p (x) . We will use p(x) and 
c(x) instead to characterize acoustic media later in full-wave simulations. Concerning the acoustic complementary 
media, we start from a layered model by transforming a single layer in acoustic wave space into two layers with 
mapping functions |1 5] . As shown in Fig.l, the original region in the x'yz Cartesian coordinate of acoustic virtual 
space with x' G [c, d] now is transformed into two regions x £ [— a, 0] and x £ [0,6] in real space xyz by mapping 
function in Fig. 1(b). For simplicity, we adopt a linear form of mapping function with /(x) : x = x' in the region 
of X G [0, b] , and g{x) : x — x' in the region of x G [—a, 0]. The boundary continuum conditions of acoustic waves 
requires /(x) and g{x) to satisfy with ^(0) = /(O) = c and g{—a) = f{b) = d at the three interfaces x = —a, x = b 
and X = 0. We obtain the constitutive parameters of the two transformed layers 1 and 2 by the following equations: 

l/pi=Ag[l/p]Al/detAg 
l/p2 = Af[l/p]Aj/detAf 

ki^kdet{Ag) 

k2^kdet{Af) ^'^^ 

WhcrcAg — diag [l/g'{x), 1, 1] — diag [—1, 1, l],Af — diag [l//'(x), 1, 1] = diag [1, 1, 1] and "det" denotes the deter- 
minant of matrix. 
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With the mapping functions satisfy to g{x) = —f{x), layer 1 and layer 2 make up of acoustic complementary media. 
And we have —pi = P2 = P and — fci = k2 = k for the mass density and bulk modulus of the two layers, where p and 
k are the constitutive parameters of the original layer. The transformation assigns exactly the same amplitude and 
phase to acoustic pressure field at the boundary x — —a and at a: = 5 as if the two layers do not exist. 

We now numerically demonstrate how to realize two-dimensional acoustic cloak with complementary media. Using a 
finite element solver (COMSOL Mutiphysics) , we calculate the acoustic scattering from a O.lm x 0.2m sized rectangle 
object incident from left to right by a monochromatic plane wave at frequency 1. 5kHz. The host media is air. The 
mass density and the sound velocity of the object are p = 2700kg/m^ and c = 5040to/s respectively, Fig. 2(a) presents 
the normalized pressure fields in colormap with respect to the incidence. In Fig. 2(b), the host media (x < 0) and 
{x > 0.2) are air, the right region (0 < a; < 0.2) is "anti-air" with parameters p = — 1.21kg /m^ and c = — 343rn/s. 
The object, embedded in the air layer to be cloaked, is exactly the same as that in Fig. 2(a). An "anti-object", 
with the same size and shape of the object at the left side, is embedded in the "anti-air" layer. The object and the 
"anti-object" are symmetrically positioned with respect to y axis. The normalized acoustic pressure fields are shown 
in Fig. 2. The regions where the normalized amplitude of the pressure fields exceeds the bound of the color bar are 
marked by white color. The absence of the scattered waves indicates invisibility of the object and the anti-object. 

To understand the perfect acoustic cloaking by complementary media in Fig. 2(b), we consider a one-dimensional 
layered system comprised of complementary media. As shown in Fig. 3, the whole structure is the combination of 
two trilayers {+A + B + A) and {—A — B — A). Layer {+A) and layer {—A) possess the same thickness dA , layer 
i+B) and layer {—B) possess the same thickness cIb- The constitutive parameters of the layers and their anti-layers 
satisfy with p+^i = — p_^,c+^ — —C-a and p+B = —p-B,c+B = —c-b respectively. According to one-dimensional 
transfer matrix method theory|,22^ a centrosymmetric trilayer {-'rA + B -\- A) can be equivalent to one homogeneous 
layer with effective characteristic impedance Reff and phase thickness <5e//- The scheme also applies for the trilayer 
{—A — B~A), its characteristic impedance -Re// ^^^^ phase thickness 5'^^^ satisfy with R'^ff = Reff and 8'^^^ = —Seff. 
The characteristic matrix of an acoustic layer has a similar form of that of transverse magnetic (TM) polarization in 
EM waves, so that we derive the characteristic matrix of the whole structure as: 



M = M1M2 = 



cos Seff iReff sin Seff 
sin Seff COS Seff 



COS Sef f 



'iReff sin Seff 
COS Sef f 
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The fact that M regresses to an identity matrix means that the complementary tri-layer as an anti-object rehabil- 
itates the whole wave field by canceling the scattering waves completely as if the whole system were nothing but the 
host media. 

We also use the concept of complementary media to demonstrate the imaging effect in acoustic waves. As a specific 
example, in Fig. 4(a) we show the pressure field distribution when a point source interacting with complementary 
media and a cylinder object. A point source with frequency 3kHz is placed at the position of (—0.2, 0) the left 
boundary of a bilayer, in which the left one is an air layer ai x G [—0.2, 0] with p — 1.21kg/m^ , c — 343m/s, and the 



right one is an "anti-air" layer at a; S [0,0.2] with p = -1.21%/? 



.3 ^_ 



-343r7i/s. A cylinder object with a radius of 



0.07m and p = 27kg /m , c — 400m/s, is located at (—0.1,0.1) in the left air layer. As shown in Fig. 4(a), we can't 
get an image of the source. If we set an "anti-object" in the "anti-air" layer right side with same size and shape of 
the cylinder object, and ensure them symmetrically positioned to axis. Then we can get an image with good quality 
as illustrated in Fig. 4(b). This implies that image effect can be achieve even an obstacle is placed in front of the 
source. 

In conclusion, we have generalized the concept of the complementary media in acoustic waves. Using this concept, 
we show examples to cloak an object outside and achieve perfect image in acoustic waves. Undoubtedly, this scheme 
can be extend to three dimensions and not constrained by the shape of the object. Recent research work on acoustic 
metamaterials have shown that bulk modulus and mass density can be designed to vary in a wide range from negative 
to positive [MJ |551 HSJ [271 HH], making transformation media attainable practically in acoustic waves as well as 
complementary media. This work was supported by the National 863 Program of China (Grant No.2006AA03Z407), 
NSFC (GrantNo.10574099, No.60674778), NECT, STCSM and Shanghai Education and Development Foundation 
(No. 06SG24). 
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FIG. 1; (a) Coordinate transformation from a single slab to bilayers. (b) The mapping function is f{x) : x — x' in the region 
of a: £ [0,6] and g{x) : x — —x in the region of a: G [— a, 0]. The transformed double layers are complementary each other when 
/(x) = —g[x) and a = b. 
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FIG. 2: Normalized acoustic pressure under an incident plane wave from left to right. (a)An object in air with mass density 
p = 27QQkg / m? and sound velocity c = 5040m/s. (b)An object in air is cloaked by complementary media. An "anti-object" 
with mass density p — —27QQkg/m? and sound velocity c — — 5040m/s is embedded in the "anti-air" region (0 < a; < 0.2). 




FIG. 3: One-dimensional acoustic layered system accounting for the perfect cloak with complementary media. (—A) and (— B) 
are the complementary media of A and B with p+A = —p~A, c+a = — c_a and p+s = —p-B, c+b ~ —c-b- 
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FIG. 4: (a) Pressure field distribution when a point source interacting with complementary media and a cylinder object. The 
cylinder object with p = 27kg/m^ ,c — 400m/s is located at the position of (—0.1,0.1). The left layer has p — 1.21kg /im? ,c = 
343m/s located at x £ [—0.2,0]. The right layer with p = —1.21kg/m?,c = — 343m/s is located at x £ [0,0.2], The host 
media {x < 0) and {x > 0.2) are air. (b) Pressure field distribution when an "anti-object" with p — —27kg/m'^,c = — 400m/s 
is embedded in the right layer at the position of (0.1,0.1). 



